Dynamical non-ergodic scaling in continuous finite-order quantum 
phase transitions 
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Abstract. - We investigate the emergence of universal dynamical scaling in quantum critical 
spin systems adiabatically driven out of equilibrium, with emphasis on quench dynamics which 
involves non-isolated critical points (i.e., critical regions) and cannot be a priori described through 
standard scaling arguments nor time-dependent perturbative approaches. Comparing to the case 
of an isolated quantum critical point, we find that non-equilibrium scaling behavior of a large 
class of physical observables may still be explained in terms of equilibrium critical exponents. 
However, the latter are in general non-trivially path-dependent, and detailed knowledge about the 
time-dependent excitation process becomes essential. In particular, we show how multiple level 
crossings within a gapless phase may completely suppress excitation depending on the control 
path. Our results typify non-ergodic scaling in continuous finite-order quantum phase transitions. 



The response of a physical system to external probes is 
an invaluable technique for unveiling the system's proper- 
ties. If the probe is dynamic, so that the Hamiltonian be- 
comes explicitly time-dependent, the system is forced out 
of equilibrium - a subject of prime practical importance 
which can soon prove full of challenges and surprises. In 
particular, understanding and manipulating the dynamics 
of zero-temperature quantum phase transitions (QPTs) [1] 
in matter has a broad significance across fields as diverse 
as quantum statistical mechanics, material science, quan- 
tum information processing, and cosmology. The extent 
to which universal quantum scaling laws persist out of 
equilibrium and encode information about the equilibrium 
phase diagram is the topic of this work. 

As early as 1970, Barouch and coworkers [2] studied the 
time-dependent T = magnetization of the anisotropic 
XY chain, and showed that equilibrium is not reached at 
the final evolution time. This non-ergodic behavior was 
later confirmed for other physical observables [3] , and the 
analysis extended to the case where the system is driven 
across its quantum critical point (QCP) by changing a 
control parameter X{t) {e.g., the magnetic field along the 
2:-axis) in time with constant quench rate r > 0. The 
emergence of non-equilibrium scaling, however, was not 
discussed. An important step was taken in Ref. [4], start- 



ing from the observation that irrespective of how slowly 
the quench occurs, adiabaticity is lost in the thermody- 
namic limit at a "freeze-out" time {tc — t) before the QCP 
is crossed. This yields a power-law prediction for the final 
density of excitations, nox(tfin) ^ ^^^{i) ^ ■, where the 
non- equilibrium critical exponent £ = dv / {vz -|- 1) is solely 
determined by the equilibrium correlation length (^) the 
dynamic critical exponents of the QCP {v and z, respec- 
tively), and the spatial dimension, d. While it is suggestive 
to realize that defect formation is a manifestation of bro- 
ken ergodicity in Barouch's sense [3], continuous experi- 
mental advances in systems ranging from ultracold atomic 
gases to quantum magnets [5] demand the applicability of 
the above Kibble-Zurek scaling (KZS) to be carefully scru- 
tinized, and the potential for more general non-ergodic 
scaling to be explored. How much information on the equi- 
librium physics is needed for reliable scaling predictions to 
be possible? 

The KZS for linear quenches across an isolated QCP 
separating two gapped phases has been confirmed by now 
for a variety of control schemes in one-dimensional (ID) 
models [6-8], including QCPs of topological nature [10] 
and noisy driving fields [11] - generalizations to repeated 
[12] and non-linear quenches [13, 14] having also been es- 
tablished. Leaving aside the case of disordered quan- 
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turn systems, where marked deviations from power-law 
behavior may be witnessed [15], the possibility of gen- 
uinely non-adiabatic scaling in low-dimensional clean sys- 
tems, whereby non-zero excitation persists (unlike KZS) 
for T — !■ oo in the thermodynamic limit, has been pointed 
out in [8]. Likewise, critical dynamics in the presence of 
non-isolated QCPs reveals a rich landscape. The need to 
modify the KZS by replacing d with the "co-dimension" 
m of the relevant critical (gapless) surface has emerged 
through a study of the 2D Kitaev model [16]. Evidence of 
non-KZS has also been reported for quenches which origi- 
nate within an extended quantum critical region [17], cross 
a multi-critical point [9, 14], or steer the system along a 
gapless critical line [18]. 

A main purpose of this work is to develop a theory and 
understanding of non-ergodic scaling for generic (power- 
law) quenches along critical regions. To achieve this goal 
two new notions are introduced, which are both path- 
dependent: one is the concept of a dominant critical point 
to establish scaling along a critical path, and the other a 
mechanism of cancellation of excitations. Besides eluci- 
dating several results recently reported in the literature, 
our analysis indicates that details on how different modes 
of excitation are accessed throughout the quench process 
are crucial. We consider several different scenarios within 
a unifying illustrative testbed, the ID anisotropic XY 
model in a transverse alternating magnetic field [19]. In 
particular, we push beyond the KZS domain - notably, by 
investigating quenches that involve a continuous Lifshitz 
QPT to a gapless phase. We also revisit the standard KZS 
and clarify how, for arbitrary continuous QPTs, it can be 
accounted for by the iterative adiabatic renormalization 
approach of Berry [20], as long as two gapped quantum 
phases are involved. Most importantly, we find that uni- 
versal dynamical scaling is obeyed by a large class of exten- 
sive physical observables throughout the quench dynamics, 
a result with practical implications in the experimental 
detection of non-ergodic scahng. 

Model Hamiltonian.— The spin-1/2 anisotropic XY 
model in a transverse alternating field is defined by [19] 
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H=- 



n- 



1-7 



(1) 

where periodic boundary conditions are assumed, that is, 
CT^ = cTq"*"^. Here, 7 S [0,1], h,5 £ [—00,00], are the 
anisotropy in the XY plane, and the uniform and alternat- 
ing magnetic field strength, respectively. This model can 
be exactly solved by following the steps outlined in [19,21]. 
The Hamiltonian (1) rewrites as H = X)/cei<'+ -^fc ~ 



keK, 



AlHkAk, where K 



+ 



{n/N, 37r/A^, ...,7r/2 



tt/N} specifies allowed momentum values, and Al. — 
{a\,a-k,b\,b-k) is a vector operator, with aj. {bj.) denot- 
ing canonical fermionic operators that create a spinless 
fermion with momentum k for even (odd) sites. Diagonal- 
ization of the reduced 4x4 Hamiltonian matrix further 
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Fig. 1: (Color on-line) Phase diagram of the alternating spin 
chain, Eq. (1). Dashed (blue) and dashed-dotted (red) lines 
define the phase boundaries for 7 — 0.5, the enclosed area cor- 
responding to the FM phase. Dashed (blue) and solid (green) 
lines correspond to 7 = 0, the enclosed area being the SF gap- 
less phase. 



yields a collection of non-interacting quasi-particles. 



H — ^ ^k,nNk,n, 



k<£K 



+ 



in terms of an appropriate number operator Nk^n for mode 
{k,n). Assuming that n labels bands in increasing energy 
order, only ek,i,ek,2 < bands are occupied at T = 0, 
with an excitation gap = e^,?, — efc,2 being given by: 



Afe(7,/i,^) 



/i^ + 52 + cos^ k + 7^ sin^ k (2) 

1/2 



2a//i2 cos2 k + 5'^{h? + 72 sin^ k) 



Quantum phase boundaries are determined by the equa- 
tions /i^ = (5^ + 1; ^2 _ _j_ rpj^g phase diagram with 

both 7 ~ 0.5 and 7 = is shown in Fig. 1. Quantum 
phases corresponding to disordered paramagnetic (PM) 
and dimer (DM) behavior emerge as depicted for arbitrary 
7. For 7 > 0, ferromagnetic order (FM phase) develops in 
the center of the phase diagram, whereas for the isotropic 
chain a superfluid (SF) phase with a gapless spectrum and 
nonbroken J7(l)-symmetry emerges. Finite-size analysis 
reveals that this model supports four distinct universality 
classes: (i) When 7 > 0, generic QCPs belong to the d = 2 
Ising universality class with critical exponents = 1, z = 
1. Different critical behavior occurs at {h ^ Q^S — ±7) 
and {h = ±1,(5 — > 0), where weak singularities in the 
ground-state energy develop (4th-order QCPs [19]), and 
v = 2, z = 1, corresponding to the alternating universality 
class; (ii) When 7 = 0, generic QCPs on the boundary 
lines belong to the Lifshitz universality class, with criti- 
cal exponents v = 1/2, z = 2. Different critical behavior 



p-2 



Dynamical non-ergodic scaling in continuous finite-order quantum phase transitions 



still occurs at {h — ±1, S —> 0), where now i/ — 1, z — 2. 
Furthermore, Ising critical exponents are recovered while 
approaching the point {h = 0,6 = 0) = O along every 
path other than {5 = 0, h ^ 0) (when there is no QCP). 
In what follows, we focus on quenching schemes where h 
and 6 are individually or simultaneously varied with time. 
We address separately different representative scenarios. 

Quenching across an isolated critical point. Sup- 
pose first that the system is linearly quenched across 
an isolated (non-multicritical) QCP that separates two 
gapped phases upon changing a single control parameter 
as SX{t) — X{t) — Xc — {t — tc)/T, where t £ [tin,ifin]- 
Without loss of generality we may assume that the sys- 
tem becomes critical at tc = 0. For finite N, the exact 
time-evolved many-body state \ip{t)) may be determined 
from numerical integration of the Schrodinger equation 
with Hamiltonian H{t), subject to |?/'(iin)} — IV^Gsl^in)), 
the latter being the ground state of H{tin). The final 
excitation density ricx(ifin) may then be computed from 
the expectation value of the instantaneous quasi-particle 
number operator over liplt)). Provided that the quench 
rate r belongs to the appropriate range^, KZS is found to 
hold irrespective of the details of the QCP and the ini- 
tial (final) quantum phase, in particular for both 2nd and 
higher-order QPTs, and independent of the path direction: 

While the excitation density is an accurate measure of 
the loss of adiabaticity in exactly-solvable models, identi- 
fying manifestations of the KZS in quantities that can be 
more directly accessible in experiments and/or meaningful 
in more general systems is essential. Remarkably, numeri- 
cal results indicate that scaling behavior holds throughout 
the quench process for a large class of physical observables, 
provided that the excess expectation value relative to the 
instantaneous ground state is considered [19]. That is, 

AO{t) = {m\0\m)~{^Gs{t)\0\ijGs{t)) 

^ ^-{^+m/i.z+i)p^(tiJ£.^^ (3) 

where /3 is a scaling exponent determined by the phys- 
ical dimension of and F is an observable-dependent 
scaling function. For instance, under a quench of the 
magnetic field strength, /i, the magnetization per site, 
Mz = obeys dynamical scaling of the 

form AM^(t) = T^~''-''''+^^/^'"'+^'^G{{t - t^) /t), whereas 
the nearest-neighbor spin correlator along the x-direction, 
XX = (X^iLi '^x'^x^^)/^^ obeys dynamical scaling of the 
form AXX{t) = r-''/('^^+i)VF((i - tc)/t), for appropriate 
scaling functions G and W , respectively - see Fig. 2. 

^That r > Tmin follows from standard adiabaticity requirements 
away from criticality, t^^^ ~ l/[™™te[ti„,tfi„]Gap{H(t))]^. The ex- 
istence of a finite upper bound Tmax follows from the fact that if 
T is arbitrarily large, a finite system never enters the impulsive 
regime, if the size-dependent contribution to the gap dominates 
over the control-dependent one. From scaling analysis under the 
assumption that the gap closes polynomially as N~'', we estimate 





0.04 














0.02 










y=8=0.5 


CO 













u ,o2 2,1/2 

h^=(8 -7 ) 














N=400 


X 


-0.02 




0.4 

0.3 








X 

< 


-0.04 


t-* 


Y=6=D.5 / 

h^iusY^ / 


% 






-0.06 


X 

- < 


0.2 




% 






-0.08 




0.1 




% 
% 





















-0.1 
-0.12 






-2 2 







-1.5 -1 -0.5 0.5 1 1.5 2 



Fig. 2: (Color on-line) Dynamical scaling under a magnetic 
field quench. Main panel: excess nearest-neighbor spin corre- 
lation per particle, AXX, vs rescaled time for the alternating 
universality class from numerical integration of the Schrodinger 
equation. Inset: excess energy per particle, AH, vs rescaled 
time for the Ising universality class from first-order adiabatic 
renormalization. 

The fact that the system becomes gapless at a single 
instant tc suggests to seek an explanation of the above 
results based on the fact that X{t) — 1/t is a small 
parameter. While a similar strategy has been imple- 
mented in [8], our emphasis is on providing a firm the- 
oretical foundation and further highlighting important as- 
sumptions. By a suitable parametrization, the relevant 
time-dependent Hamiltonian may be written as H{t) = 
He + [A(t) - Xc]Hi = Hc + {t- tc) /tHi , with He quantum- 
critical in the thermodynamic limit. Let {Em{t)} and 
{|V'm(i))} denote the snapshot eigenvalues and (orthonor- 
mal) eigenvectors of H{t), where |V'o(i)) = \^Gs{t)) and 
assume that: (i) no level crossing occurs throughout the 
evolution; (ii) the derivatives of all the spectral projectors 
{\'ipm{t)){'4'm{t)\} are sufficiently smooth. The (normal- 
ized) time-evolved state reads 

|V(t)) - c^{t)\Mt)) + c™(i)IV'm(i)), 

for coefficients to be determined. Since for a truly adi- 
abatic evolution no excitation is induced in spite of the 
fact that the eigenstates of H{t) evolve in time, ap- 
propriately subtracting (following Berry, 'renormalizing') 
the adiabatic contribution is essential for quantifying 
the leading non-adiabatic correction. This is achieved 
in two steps [22]: (i) Effect a canonical transformation 
to a 'comoving frame,' where in the zeroth-order adia- 
batic limit T — > oo the comoving state vector |V'(0) = 
U{t;ti-a)\%lj{ti-a)) is frozen up to a phase factor, that is, 
■!/'(i)) = e^*^"'*) IV'o(^in)), where To{t) includes in general 
both the Berry phase and the dynamical phase; (ii) Evalu- 
ate the first-order correction to the comoving-frame prop- 
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agator via Dyson series expansion. Transforming back 

to the physical frame, Cm{t) = (V''m(iin)|C^(i; iin)IV'o(im)), 
to first-order in A we finally obtain (in units h = 1), 



-iTo(t) 



-hO(A2), and 



.(1) 



it) 



-iTUt) 



JUn 



Emit') - Eoit') 



Am(t) = Emit) - Eoit) . 



(4) 



Knowledge of the time-dependent state enables arbi- 
trary physical quantities of interest to be computed, in 
particular the total timc-dcpcndcnt excitation probabil- 
ity Pex(t) = Em^oMt)\\ Given Eq. (4), the latter 
formally recovers the expression given in [8], which cap- 
tures the contribution to the density of excitations from 
states directly connected to \tpo{t)) via Hi ^. Dynami- 
cal scaling emerges once the above result is supplemented 
by scaling assiimptions on three fundamental dynamical 
variables: the time-dependent excitation energy above the 
ground state; the time-dependent matrix elements of the 
perturbation; and the density of excited states, p{E), at 
the energy scale characterizing adiabaticity- breaking, 
which allows to change discrete sums over excited states 
to integrals. That is, close to the QCP we assume that: 

Emit) - Eoit) = 5\it)'''fmi/^mito)/5\itr'), 

{^mit)\Hi\Mt)) = 6Xit)-^-'gmiAmitc)/SXit)n, 

piE) ~ E'''^-\ (5) 

where the scaling functions fm,9m satisfy i) igm) is 
constant when x ^ 0; ii) (<?„) c< x when x ^ oo 
^. Having the scaling assumptions at hand, integration 
over excited states is performed by moving to dimension- 
less variables ( = [t - tc)/i = {t - tc)r-''^/(''^+i) and 
ri = A,„(tc)t = A„(ic)r''^/('^^+i). Since at the QCP 
the integrand in Eq. (4) develops a simple pole, while the 

phase g^^fin''''^'"^''' becomes stationary, contributions away 
from the QCP may be neglected, allowing the desired scal- 
ing factor to be isolated, up to a regular function depend- 
ing only on C^. Thus, the scaling of the excitation density 
and diagonal observables such as the residual energy is di- 
rectly determined as nex(C) = t-'*''/('^^+i)S(C), Aif(C) = 
^-{d+z)i'/(uz+i)'Y(^(^'j^ g^jgQ pjg_ 2. For a generic observ- 
able, if the additional scaling condition 

{Mt)\ Hmit)) = 5Xit)^qmiAmitc)/SXit)n, (6) 

holds for all the excitations m involved in the process 
for an appropriate scaling function q„i, then AO ~ 
^-{ud+f3)/{uz+i) _ consistent with Eq. (3). 



^In particular, since one-body perturbations Hi are considered in 
the present analysis, the first-order excitation probability, F^]} (t) = 
Yl^rnj^o {'^■m (t)\^ , Coincides with the single-mode quasiparticle contri- 
bution, {Ni^ „) = 1, to the total time-dependent excitation density. 

^Notc that p ~ (-'^/E, with ^'^ ~ ^m^d-m ^ ^-m/zi^d-m foj. ^ 

{d — m)-dimensional critical surface. 



Two remarks are in order. First, the above argument 
directly explains the dynamical scaling reportc;d in [19] 
for generalized entanglement relative to the fermionic al- 
gebra u(iV) [23], whose ground-state equilibrium behav- 
ior directly reflects the fluctuations of the total number 
operator. Second, the derivation naturally extends to 
a generic non-linear power-law quench, that is, 5Xit) = 
A(t) - Ac |(i - tc)/T|"sign(i - t^), a > 0. Provided that 
the typical time scale for adiabaticity breaking is rede- 
fined as ia ~ T°"^^/^^~^°'"^\ the same scaling assumptions 
in Eqs. (5)- (6) lead to dynamical scaling behavior of the 
form Hex ~ r-"''''/(«''^+i), and AO ~ T--oc{d,^+fi)/ia'^z+i) ^ 
throughout the whole time evolution 

Quenching along paths involving a finite number 
of critical modes.— A first situation which is beyond the 
standard KZS discussed thus far arises in quenches that 
force the system along a critical line, yet are dominated 
by a finite number of participating excitations. Formally, 
this makes it possible to obtain the non-equilibrium ex- 
ponent for riex through application of the KZS, provided 
care is taken in defining the static exponents through a 
limiting path-dependent process where, along the quench 
of interest, a simultaneous expansion with respect to both 
the control parameter and the relevant critical mode(s) 
is taken. Consider a quenching scheme where both h 
and 6 arc c;hangcd according to t/r while 7 = (path 
F ^ O -* G in Fig. 1). While Eq. (2) shows that 
the mode k = 7r/2 is critical throughout the process 
(A^/2(i) = 0, for all t e [tin,ifin] as N ^ oo), numer- 
ical data indicate that excitation sets in only when the 
point O is passed, see Fig. 3. As remarked, the static crit- 
ical exponents at O are z = = 1, which differ from 
the critical exponents (2; = 2, = 1/2) of all other critical 
points along this line. Indeed, the non-equilibrium expo- 
nent is solely determined by the static exponents of this 
QCP along the chosen path, riex ~ r-^'/^'^^+i) = t-^/^. 
We term a QCP which belongs to a different universal- 
ity class than all other critical points along a critical line 
and sets the non-ergodic scaling a dominant critical point 
for that line ®. Physically, although A^/2 closes along 
the critical line in the thermodynamic limit, a level cross- 
ing which brings all bands together only occurs at O - 
still allowing the time-evolved state to adiabatically fol- 
low the snapshot ground state until then. The follow- 
ing independent confirmations may be invoked in support 
of the above argument. First, consider the anisotropic 
quench process analyzed in [18], whereby j{t) is changed 
linearly along the critical line = 5^-1-1. By Taylor- 
expanding Afe in Eq. (2) around fc = 0, 7 = reveals 
that v = l,z ~ 2 a.t the dominant QCP (7 = 0,h,5), 
whereas i^=l,z = lfor77^0 along the line. Accord- 



*The perturbative derivation as presented strictly applies to 
quenches across an isolated QCP which is not multi-critical. We 
defer application of the perturbative derivation to a multi-critical 

point to a forthcoming analysis. 

'The point O is multi-critical. However, quenches across a multi- 
critical point need not satisfy KZS, see [14]. 
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Fig. 3: (Color on-line) Main panel: dynamical scaling of the 
excitation density for a simultaneous linear quench of h and S 
along the gapless critical line F O ^ G. Inset: log-log plot 
of the final excitation density vs r along the path A' D ^ O. 

ingly, riox ^ t~^^'^. While this coincides with the result 
obtained in [18], the underlying physical explanation is dif- 
ferent. Plots of the rescaled excitation density (ric^T^^^) vs 
the rescaled time (t/r^/^) would collapse onto one another 
for different r within the appropriate range, in complete 
analogy with Fig. 3. Second, loss of adiabaticity at a sin- 
gle point can also explain the scaling behavior observed 
for a AFM-to-FM quench (or a critical-to-FM quench) in 
the XXZ model [17], whereby the control path involves the 
gapless critical region — 1 < A < 1 and the dominant crit- 
ical point A = 1 belongs to a different universality class. 
Lastly, the concept of a dominant QCP remains useful for 
a power-law quench, which leads to the scaling behavior 
nex T-adv/iai^z+i) ^ with v and z being the critical ex- 
ponents of the dominant QCP along the critical line. 

Quenching along paths involving an infinite 
number of critical modes.— More complex scenarios 
emerge when uncountably many modes of excitations can 
compete during the quench. Focusing on the isotropic 
limit 7 = 0, we contrast two representative situations 
where the Lifshitz QPT is involved: (I) Magnetic quenches 
along the path D ^ E (PM SF ^ PM); (II) Al- 
ternating quenches along the path A ^ B ^ C (DM — > 
SF DM). Since [M^,H] = 0, in both cases the allowed 
excitation must comply with a non-trivial dynamical con- 
straint. Along the path D ^ O ^ E, this forces the final 
state to be the same as the initial ground state up to a 
global phase factor, leading to nex(ifin) ~ Although for 
a magnetic quench this may be viewed as a consequence 
of the fact that the dynamics simply acts as a relabelling 
of the snapshot eigenstates, the same scaling holds for any 
quench which begins or ends in the gapless phase - for in- 
stance, a (5-quench along the path A ^ B. Because these 
quenches take the system through a critical line in mo- 
mentum space, d — m = 1 (as opposed to d — to = for 



an isolated QCP), the observed scaling is consistent with 
the recent prediction nox(tfin) ^ t-""'/(''^+i) [16]. 

One may naively expect the same scaling to hold for 
path (II), which also connects two gapped phases, albeit 
different than in (I). Unlike in the standard KZS, how- 
ever, details about the initial and final phases as well as 
the time-dependent excitation pattern become important. 
Specifically, along path (II) we find nox(ifin) t^^/^. An 
explanation may be obtained by exploiting the fact that 
due to [/(l)-symmetry, the fermion number is conserved. 
This allows the reduced 4x4 matrix Hk to be decoupled 
into two 2x2 matrices by interchanging the order of the 
basis vectors a^k and fej,. Thus, H±k = Wl.f.H' ±kW±k, 
where wl = ial,bl), W^^ = {a^k,b^k), and 

). ,7, 

For such a two-level system, the asymptotic excitation 
probability may be computed from the Landau-Zener 
transition formula [24], yielding pk ~ g-27rcos kr ^ Upon 
integrating over all modes, we find 

1 /""/^ 

ncx(ifin) = - / dkpk ~ (8) 

Note that because pk is independent on h, the result in 
Eq. (8) may be interpreted as implying that traversing 
the gapless phase produces the same excitation density 
as crossing the single QCP O by translating path (II) at 
h = 0, which determines the non-equilibrium exponent. 

Physical insight into what may be responsible for the 
different behavior observed in the two quenches is gained 
by looking at the excitation spectrum along the two paths. 
Notice that, once the energy eigenvalues are specified at 
the initial time, (efc,i(iin) < efc,2(iin) < efe,3(iin) < (-kAiUn) 
in our case), the same relative ordering need not hold at 
the final time if a level crossing is encountered during the 
quench - see Fig. 4(a), (b). In the critical region, a pair 
of modes {k,n) and {k,n') undergoes a level crossing if 
h? — 5^ ^ cos^ k. If the number of such level crossings for 
fixed n, n' is even, the net contribution to the final exci- 
tation from momentum k is zero, since the final occupied 
bands are the same as the initial ones - see Fig. 4(c). No 
cancellation is in place if either an odd number of level 
crossings from the same pair or if different pairs (n, n') 
are involved. The latter situation is realized for all k along 
path (I) (/i-quench. Fig. 4(a)) and also for the path A ^ B 
((5-quench). For a 5-quench along path (II), the net exci- 
tation from the gapless phase turns out to be completely 
canceled (as seen in Fig. 4(d), where the quench starts and 
ends symmetrically within the gapless phase). This only 
leaves the two boundary critical lines h = ±5 as contribut- 
ing to the excitation, thus a finite set of critical modes (a 
single one in fact, k — tt/2, see Fig. 4(b)). Interestingly, a 
similar cancellation mechanism was verified for repeated 
quenches across an isolated QCP [12]. While a thorough 
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Fig. 4: (Color on-line) Panels (a — c): band structure for 
different momentum modes k vs quench parameter {i.e., h or 
5). Band ordering is determined by band index n — 1,2,3,4 
at initial time t = tin, whereas dash-dotted (black) < dotted 
(blue) < dashed (green) < solid (red) at a generic time t along 
the control path. Panel (d): final excitation density, ncx(tfln), 
vs 5 for a quench of the alternation strength 5 at fixed h = 0.5 
within the critical region. For all cases, TV = 400. 



analysis is beyond our current scope, we suggest that even 
participation from the same (pair of) snapshot excitations 
may be at the root of this cancellation in both scenar- 
ios. Here, we further test this conjecture by examining 
the path A' ^ D ^ O, for which an effective two-level 
LZ mapping is no longer possible. Unlike A ^ B ^ C, 
two intermediate phases are now crossed, and the initial 
and final phase differ from one another, yet analysis of 
^k,n(t) reveals that the two paths are equivalent in terms 
of participation of critical modes. Numerical results con- 
firm that 7icx(ifin) ~ T~^/'^, Fig. 3. 

Conclusions.— Non-ergodic dynamical scaling is fully 
captured by first-order adiabatic renormalization for suf- 
ficiently slow quenches involving a simple isolated QCP. 
Beyond this regime, we find that non-equilibrium expo- 
nents remain expressible by combinations of equilibrium 
(path-dependent) ones in all the scenarios under examina- 
tion, however a detailed characterization of both the static 
phase diagram and the accessible low-lying excitations is 
necessary for quantitative predictions. Ultimately, scaling 
behavior appears to be the same for control paths which 
share an equivalent excitation structure. While yet differ- 
ent non-ergodic scaling may arise in more complex systems 
{e.g., infinite-order Berezinskii-Kosterlitz-Thouless QPTs 
[17] as well as models with infinite coordination [25]), a 
deeper analysis of how competing many-body excitations 
contribute and interfere during a quench may shed further 
light on non-equilibrium quantum critical physics. 
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